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Abstract  Hlgh-resolution  global  ocean  models  forced  by  atmospheric  fields  and  tides  are  beginning  to 
display  realistic  internal  gravity  wave  spectra,  especially  as  model  resolution  increases.  This  paper  examines 
internal  waves  in  global  simulations  with  0.08°  and  0.04°  (~8  and  4  km)  horizontal  resolutions,  respectively. 
Freq  uency  spectra  of  internal  wave  horizontal  kinetic  energy  in  the  North  Pacific  lie  closer  to  observations  in  the 
0.04°  simulation  than  in  the  0.08°  simulation.  The  horizontal  wave  number  arx)  frequency  [K-co)  kinetic  energy 
spectra  contain  peaks  in  the  semidiurnal  tidal  band  and  near-inertial  band,  along  with  a  broadband  frequency 
continuum  aligned  along  the  linear  dispersion  relations  of  low-vertical-mode  internal  waves.  Spectral  kinetic 
energy  transfers  desaibe  the  rate  at  which  nonlinear  mechanisms  remove  or  supply  kinetic  energy  in  specific 
K-(o  ranges.  Energy  is  transferred  out  of  low-mode  inertial  and  semidiurnal  internal  waves  into  a  broad 
continuum  of  higher-frequency  and  higher-wave  number  internal  waves. 


1.  Introduction 

In  this  paper,  we  examine  the  extent  to  which  high-resolution  global  ocean  models  that  are  simultaneously 
forced  by  atmospheric  fields  and  tides  [e.g.,  Arbicetal.,  2010, 2012a;  Muller  et  al.,  2012]  are  able  to  simulate 
the  internal  gravity  wave  spectrum  and  the  extent  to  which  nonlinear  internal  wave-wave  interactions 
contribute  to  the  simulated  sprectrum.  For  about  20 years,  basin-  and  global-scale  ocean  models  have  been 
able  to  simulate  mesoscale  eddies  with  some  degree  of  realism  [McClean  et  aL,  1997].  The  spectrum  of 
internal  gravity  waves  (hereafter  referred  to  simply  as  Internal  waves")  represents  a  new  frontier  for  global 
models.  Because  internal  waves  redistribute  large  amounts  of  energy  over  long  distances,  analyses  of  the 
global  distribution  of  internal  waves  and  their  interactions  are  relevant  to  our  understanding  and 
simulation  of  large-scale  ocean  dynamics  [Munk  and  Wunsch,  1998]. 

Internal  wave  space-time  characteristics  in  deep  water  have  been  represented  by  a  nearly  universal  spectrum 
[Garrett  and  Munk,  1975;  Cairns  and  Williams,  1976].  The  internal  wave  spectrum  covers  a  range  of  scales — in 
time  from  the  buoyancy  period  (~10min  in  the  thermocline)  to  the  inertial  period  (~1  day  but  depends  on 
latitude),  in  vertical  scale  from  ~1  m  to  the  ocean  depth,  and  in  horizontal  scale  from  ~  100  m  to  greater 
than  100km.  The  internal  wave  spectrum  is  fed  primarily  by  near-inertial  waves  arising  from  surface  wind 
forcing  [e.g,  D'Asaro,  1984;  S'dverthorne  and  Toole,  2009;  Simmons  and  Alford,  2012]  and  by  internal  tides 
arising  from  barotropic  tidal  flow  over  topography  [e.g,  Garrett  and  Kunze,  2(X)7].  Nonlinear  internal  wave- 
wave  interactions  have  been  proprosed  as  the  primary  mechanism  to  transfer  energy  out  of  the  generation 
frequencies  and  wave  numbers  to  the  broader  range  of  frequencies  and  wave  numbers  in  the  internal 
wave  sprectrum  [Muller  et  al.,  1986;  Polzin,  2004].  In  an  idealized  two-dimensional  model,  with  sjsatial 
dimensions  in  the  horizontal  and  vertical  that  was  forced  at  the  near-inertial  and  semidiurnal  frequencies, 
Sugiyama  et  al.  [2009]  found  that  internal  wave-wave  interactions  filled  out  a  frequency  wave  number 
sp)ectrum.  In  the  present  study,  we  examine  the  internal  waves  and  their  interactions  in  realistic  three- 
dimensional  global  simulations. 

A  review  of  nonlinear  interactions  in  the  ocean  internal  wave  field  by  Muller  et  al.  [1986]  focuses  on  early 
theoretical  work.  The  more  recent  review  by  Polzin  et  al.  [2014]  focuses  on  the  development  and 
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application  of  a  practical  fine-scale  parameterization  for  turbulent  dissipration  and  mixing  based  on  the 
energy  cascade  through  the  vertical  wave  number  sprectrum  to  small  vertical  scales.  This  fine-scale 
jsarameterization  has  found  wide  use  observationaliy. 

Instead  of  isolating  sp)ecific  triads,  as  in,  for  instance.  Sun  and  Pinkel  [2013]  and  MacKinnon  et  al.  [2013], 
we  will  calculate  the  change  in  the  horizontal  kinetic  energy  sp)ectrum  due  to  nonlinear  momentum 
advection,  thus  capturing  the  full  range  of  nonlinear  effects  at  the  exprense  of  identifying  sprecific 
interactions  in  detail.  We  will  compute  horizontal  kinetic  energy  spjectra  E  and  nonlinear  sprectral 
kinetic  energy  transfers  T  in  the  horizontal  wave  number-frequency  (K-to)  domain  based  on  surface 
velocities  of  realistic  high-resolution  global  ocean  circulation  models  forced  with  atmospheric  fields  and 
tides.  Note  that  K=(k^  +  /^)'^^  is  the  horizontal  wave  number,  with  k  and  /  being  the  zonal  and 
meridional  wave  numbers,  resprectively.  As  discussed  further  below,  our  focus  on  horizontal  wave 
number  is  due  to  the  availability  of  surface  model  fields  and  the  computational  difficulty  of  prerforming 
sprectral  computations  with  fully  three-dimensional  high-resolution  model  output.  Our  K-a>  internal  wave 
sprectral  transfers  are  akin  to  transfers  that  have  long  been  computed  In  the  horizontal  wave  number 
domain  for  quasi-geostrophic  (QG)  motions  [e.g,  Scott  and  Wang,  2005].  Arbic  et  al.  [2012b,  2014] 
extended  the  sprectral  analysis  of  QG  kinetic  energy  transfers  to  the  frequency-  and  wave  number- 
frequency  domains. 

In  section  2,  we  introduce  the  model  output  and  observational  data  used  in  this  paper.  In  section  3,  the 
derivation  of  kinetic  energy  spectra,  £,  and  sjsectral  energy  transfers,  T,  in  [K-oi]  spiace  is  presented,  in 
section  4,  we  compare  the  simulated  kinetic  energy  frequency  spjectra  with  spectra  from  current  meter 
observations.  In  section  5,  we  demonstrate  that  the  modeled  E  and  T  in  (K-ca)  spjace  are  concentrated 
along  the  dispjersion  relations  of  low  vertical  modes.  In  section  6,  we  interpret  the  (K-w)  spjectra  in  terms  of 
triad  interactions.  A  summary  is  presented  in  section  7. 

2.  Model  Output  and  Observational  Database 

We  utilize  output  in  a  North  Pacific  domain  spjanning  29°N-43°N  and  150°W-176°W  from  two  global 
simulations  of  the  HYbrid  Coordinate  Ocean  Model  (HYCOM)  [Chassignet  et  al.,  2009].  The  model  is 
hydrostatic  with  32  layers  in  the  vertical  direction.  The  chosen  analysis  domain  is  north  of  the  Hawaiian 
islands  in  a  region  of  substantial  internal  tide  activity.  The  nominal  horizontal  resolutions  of  the 
simulations  are  0.08°  (HYCOM12)  and  0.04°  (HYCOM25),  correspjonding  to  horizontal  grid  spjacing  of  ~8 
and  4km,  respectively.  Note  that  the  vertical  resolution  is  the  same  for  both  model  simulations.  The 
models  are  forced  by  astronomical  tides  and  3-hourly  atmospheric  fields.  As  discussed  in  Arbic  etal.  [2010], 
the  models  employ  a  pjarameterized  topjographic  internal  wave  drag  and  a  scalar  approximation  for  the 
self-attraction  and  loadirrg  term.  The  spectral  calculations  in  this  pjapjer  use  6  (3)  months  of  hourly  output 
from  HYCOM12  (HYCOM25). 

The  internal  wave  spjectrum  is  typically  described  in  terms  of  frequency  and  vertical  modes  [e.g,  Garrett  and 
Munk,  1975].  However,  the  volume  of  three-dimensional  output  from  the  HYCOM12  and  HYCOM25 
simulations  is  extremely  large.  Therefore,  for  computational  feasibility,  we  have  chosen  to  restrict  this  first 
analysis  of  HYCOM  internal  wave  spjectra  to  surface  horizontal  velocities  rather  than  to  three-dimensional 
velocity  fields.  The  dispjersion  relations  of  low  vertical  modes  will  be  identified  in  the  horizontal  wave 
number-frequency  spjectra. 

The  lower-resolution  HYCOM12  has  been  extensively  compjared  to  available  current  meter  measurements  of 
the  global  three-dimensional  tidal  velocity  field  [Timko  et  al.,  2013]  and  to  satellite  altimetiy<onstrained 
estimates  of  the  barotropic  and  baroclinic  compjonents  of  surface  tidal  elevations  [Shriver  et  al.,  2012].  The 
latter  study  computed  a  7.0  cm  global  root-mean-square  M2  sea  surface  elevation  error  in  the  deep  ocean 
with  respect  to  the  TPX08-atlas  [Egbert  et  al.,  1 994].  This  error  is  compjarable  to  errors  in  other  forward  tide 
models  employing  a  parameterized  topjographic  internal  wave  drag.  The  frequency  spjectra  of  kinetic 
energy  f(to)  from  the  model  are  compared  to  frequency  spjectra  from  historical  velocity  observations  at 
seven  moorings  within  our  North  Pacific  domain  along  175°W  and  152°W.  The  moorings  were  deployed  as 
pjart  of  the  Pacific  Zonal  Exploration  Array  [Schmitz,  1988].  The  instruments  are  located  in  the  uppjer  ocean 
between  143  and  168  m  with  record  lengths  varying  between  348  and  705  days  and  sampling  pjeriods 
ranging  from  1 5  min  to  1  h. 
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Figure  2. The  Km  spectra  of  kinetic  energy  E(K,  m)  {(m/s)^[(d) 
(km)]}  on  a  logio  scale  for  (a)  HYCOM12  and  (b)  HYCOM25.  The 
white  curves  represent  bounding  linear  dispersion  relations  of 
the  first  three  internal  \Arave  vertical  modes  (solid  =  mode  1 ; 
dashed  =  mode  2;  and  dash  dotted  =  mode  3).  The  mode  bands 
are  labeled  inside  the  panel  with  white  lettering.  Below  and  to  the 
right  of  the  color  contoured  E{K,  m),  the  spectra  E(m)  integrated 
over  all  horizontal  wave  numbers  K,  and  EfnviK)  integrated  over 
frequencies  larger  than  the  inertial  frequency  (q>>  /jqn)-  are 
shown,  respectively.  In  the  £(03)  spectra,  cyan  lines  indicate  the 
inertial  frequencies  at  the  bounding  latitudes  29°N  and  43°N, 
while  M2  and  M4  frequencies  are  also  indicated. 


3.  Derivation  of  Kinetic  Energy 
Spectra  and  Nonlinear 
Spectral  Transfers 

We  compute  E  and  T  directly  from  the 
surface  velocity  without  employing  the  (JG 
assumption  used  in  Scott  and  Wang  [2(X)5] 
and  Arbic  et  at.  [2012b,  2014].  We  begin  with 
the  momentum  equation 

^+[(uV)u]  =  OT,  (1) 

where  OT  refers  to  other  terms,  not  considered 
in  the  present  analysis,  and  u  =  (u,  v)  is  the 
horizontal  velocity  vector  with  zonal  and 
meridional  components  u  and  v,  respiectively. 
Note  that  in  equation  (1),  we  exclude  nonlinear 
terms  of  the  fomn  wdu/dz,  where  w  is  the 
vertical  velocity.  Taking  the  Fourier  transform 
of  equation  (1),  and  multiplying  by  the 
complex  conjugate  of  the  Fourier  transform  of 
the  horizontal  velocity  vector,  yields  a  spiectral 
energy  equation 

|£(K,  m)  =  T{K,  0))  +  OT.  (2) 

Thus,  £  and  T are  defined  by 

E(K,m)=  (3) 

T(K,  ra)  = /?e{  l/••[u  •Vu]},  (4) 

where  *  denotes  complex  conjugate,  *  refers  to 
the  Fourier  transform,  and  Re  refers  to  the  real 
part  of  a  complex  number.  From  equation  (2), 
negative  (pxisitive)  values  of  T[K,  co)  imply  that 
nonlinear  interactions  represent  a  sink  (source) 
of  kinetic  energy  for  the  given  (K,  to).  Before  the 
sfiectral  analysis  is  performed,  the  time  series 
are  detrended  and  windowed  in  space  and 
time  following  standard  procedures  of  time 
series  analysis. 

4.  Comparison  With  Observations 

Near-inertial  p>eaks  in  the  £(to)  kinetic  energy 
spectra  at  the  seven  mooring  locations  are 
reasonably  consistent  between  the  models 


Figure  1.  (a  g)  Frequency  spectra  flo))  of  surface  kinetic  energy  {(m/s)^(d)}  from  moored  current  meter  observations  (red), 
HYCOM12  (blue),  and  HYCOM25  (green)  and  the  GM76  [Cairns  and  Williams,  1976]  model  spectrum  (dashed  black).  The 
spectra  are  logarithmically  smoothed  for  display  purposes.  GM76  spectra  are  computed  with  reference  values  of  total 
energy  of  the  internal  wavefield  and  stratification  set  to  Eg  =  63e  ^m^s  ^andNo  =  53e  ^s  ',  respectively.  Vertical  lines 
are  drawn  at  f,  2f,  and  f+M2,  while  dashed  vertical  lines  are  drawn  at  the  four  semidiurnal  (M2,  S2,  N2,  and  K2)  and  four 
diurnal  [K^,  O1,  Pi,  and  Q,)  tidal  frequencies  forced  in  the  model,  (h)  The  analysis  region  in  the  northeast  Pacific  The  blue 
stars  denote  the  locations  of  moored  observations  used  in  Figures  la  1g. 
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Figure  3.  As  in  Figure  2  but  for  nonlinear  spectral  kinetic  energy 
transfer  r(ff,  ro)  {(10“®  WAg)  [(d)  (km)]},  with  linear  colorbar  and 
dispersion  curves  artd  mode  band  numbers  shown  in  black. 
Energy  is  taken  out  of  the  blue  regions  (low  mode  near  inertial 
and  semidiurnal  internal  waves)  and  added  to  the  red  regions 
(higher  frequency  and  higher  wave  number  internal  waves). 


and  observations  (Figures  la-lg)  except  for  the 
observation  at  41.94°N,  151.92°W,  where  the 
HYCOM25  and  HYCOM12  peaks  are  3  and  9 
times  larger,  respectively,  than  the  observed 
peak  (Figure  1e).  This  is  likely  because  of 
different  winter  storm  forcing  during  the  year 
of  the  current  meter  measurements  than  in 
the  model  [D'Asaro,  1985].  In  the  other 
locations,  the  near-inertial  peaks  differ  by 
factors  ranging  between  0.3  and  2.5  and  thus 
are  within  an  order  of  magnitude  of  each 
other.  Oiumal  peaks  are  visible  in  the 
observations  at  the  eastern  but  not  western 
location  and  not  in  the  models  at  any  location. 
Semidiurnal  tidal  peaks  are  visible  in  all 
observational  and  model  spectra  shown  in 
Figure  1.  Above  the  semidiurnal  tidal 
frequency,  the  energy  in  the  observations 
decreases  approximately  as  consistent 
with  the  GM76  [Cairns  and  Williams,  1976] 
spectra.  The  models,  especially  HYCOM12, 
have  a  steeper  roll-off  at  high  frequencies.  At 
super-tidal  frequencies,  I-IYCOM25  rolls  off  less 
steeply,  and  thus  lies  closer  to  the 
observations  than  HYCOM12.  We  speculate 
that  this  is  due  to  better  resolution  of  higher 
vertical  modes,  which  are  required  to  fill  out 
the  internal  wave  spectrum  [Garrett  and  Munk, 
1975].  At  some  locations,  HYCOM25  displays 
peaks  at  frequencies  of  2f  and  f+  eoui  (f  refers 
to  the  Coriolis  frequency  and  (0^2  refers  to  the 
M2  tidal  frequency)  that  are  also  visible  in 
the  observations. 

5.  Horizontal  Wave  Number- 
Frequency  Spectra  and 
Spectral  Transfers 


The  spectrum  E[K,  co),  along  with  the  dispersion 
cun/es  of  internal  waves,  are  shown  in  Figure  2. 
To  compute  the  disprersion  curves,  we  solve  for 
the  eigenspeeds  €„  of  the  nth  vertical  mode 
[Munk,  1981]  using  a  shooting  method  code  (Glenn  Flierl,  personal  communication,  1995)  applied  to  the 
internal  wave  Sturm-Liouville  problem.  Stratification  profiles  N^{z),  computed  from  time-averaged 
temperatures  and  salinities  at  individual  grid  points  from  the  I-IYCOM25  output,  are  taken  as  inputs  to 
the  Sturm-Liouville  problem.  The  dispersion  relation  for  a  prarticular  vertical  mode  is  given  by 

=f^  +ciK^.  (5) 


For  the  analysis  region,  f  varies  with  latitude  and  W^(z)  varies  from  one  model  grid  point  to  the  next  We 
computed  eigenspeeds  c„  for  each  of  the  first  three  vertical  modes  at  all  model  grid  points  along  the  29°N 
and  43°N  bounding  latitudes.  The  bounding  dispersion  curves  for  each  mode  are  computed  using  the 
minimum  Cn  value  along  the  northern  latitude  and  the  maximum  c„  value  along  the  southern  latitude. 
Because  all  computations  made  from  these  high-resolution  models  are  time  consuming,  we  did  not 
separately  compute  the  dispersion  curves  for  the  HYCOM 12  stratification  but  use  those  from  HYCOM25. 
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The  models'  K-ca  sp)ectra  (Figure  2)  display  lobes 
between  the  bounding  dispersion  curves. 
Because  of  the  higher  resolution  of  HYCOM25, 
its  kinetic  energy  spectrum  E{K,  m)  has  higher 
energy  levels  extending  to  smaller  spatial  and 
temporal  scales  than  HVCOM12.  Both  models 
exhibit  the  largest  spectral  kinetic  energy  in 
the  inertial  frequency  range,  defined  by  the 
29’N-43°N  bounding  latitudes  of  the  analysis 
region.  Further,  energy  accumulation  is  found 
along  the  dispersion  relations  for  higher  vertical 
modes.  The  peaks  for  high  vertical  modes  are 
centered  at  horizorrtal  wavelengths  of  about 
50  km  (0.1 25  rad/km)  atxl  37  km  (0.1 70  rad/km) 
for  FIYCOM12  and  HYCOM25,  respectively.  In 
the  E{K,  a>)  spectra,  the  semidiurnal  tidal  lines 
are  visible  at  around  12rad/d  and  peak  within 
the  bcxinding  dispersion  relations  for  mcxies 
1,  2,  and  3  with  approximate  wavelengths  of 
144,  72,  and  50  km,  respectively.  However, 
substantial  energy  is  found  between  those 
peaks,  not  covered  by  the  range  of  mcxial 
dispersion  relations.  The  spectra  £(cu)  integrated 
over  all  wavelengths  also  display  semidiurnal 
peaks.  Tidal  harmonics,  such  as  M4.  are  seen  in 
the  HYCOM12  spectrum.  Diurnal  tides  are  not 
readily  apparent  in  the  frequency  spectra  E{(o). 
In  ackiition,  we  present  kinetic  energy  spectra 
frrwffO  integrated  over  all  frequencies  higher 
than  the  inertial  frequency  at  the  lower 
boundary  of  the  region.  Thus,  ftrw(fO  represents 
the  horizontal  wave  number  spectra  for  internal 
wave  horizontal  kinetic  energy.  In  EmM,  a 
peak  is  visible  at  the  0.044  rad/km  wave  number 
corresponding  to  the  rrxxie  1  internal  tide. 

The  nonlinear  spectral  kinetic  energy  transfer 
UK,  m)  is  displayed  in  Figure  3.  The  largest 
negative  (sink)  values  of  the  nonlinear 
spectral  kinetic  energy  transfer  T{K,  or)  show  energy  being  transferred  out  of  the  low-wave  number 
near-inertial  and  supermode  3  semidiurnal  bands  in  HYCOM25  (Figure  3b),  the  latter  coinciding  with 
the  large  high-mode  E{K,m)  values  discussed  earlier.  The  negative  T{K,m)  values  imply  that 
nonlinearities  extract  energy  from  these  regions  of  the  K-co  spectrum  to  feed  other  spectral  regions. 
The  largest  positive  (source)  TIK.  m)  shows  energy  being  transferred  into  the  supermode  3  near-inertial  arxf 
mode  2  2f  bands  in  HYCOM25.  As  with  the  spectra  E(K,  oj),  large  UK,  eo)  extends  to  higher  wave  numbers  and 
frequencies  in  HYCOM25  than  HYCOM12,  demonstrating  that  nonlinear  interactions  are  more  active  in 
HYCOM25.  Also  as  in  E(K,ca),  the  semidiurnal  internal  tides  are  visible  as  lines  in  TiK.ai)  and  H<a), 
with  segments  of  increased  energy  transfer  coinciding  with  the  ranges  given  by  the  respective  dispersion 
relations. 

Along  the  29°N  southern  boundary  of  the  domaia  the  internal  wave  transfer  rates  Tnw(fO  has  a  large 
negative  (sink)  peak  at  wavelengths  larger  than  approximately  300  km.  The  dominant  motions 
contributing  to  this  sink  are  near  inertial.  For  wavelengths  less  than  300  km,  Tnw(fO  is  characterized  by  a 
small  positive  (source)  peak  followed  by  a  larger  negative  (sink)  peak  at  still  higher  wave  numbers.  While 
the  peak  at  wavelengths  larger  than  300  km  is  consistent  between  HYCOM12  and  HYCOM25,  the 


a 


Figure  4.  Sketch  of  internal  wave  wave  triad  interactions  in  /Cq> 
space,  (a)  The  dispersion  relations  of  four  internal  wave  vertical 
modes.  The  cross  section  of  the  green  dashed  and  yellow  lines 
illustrates  the  triad  interactioa  where  mode  1  arxJ  mode  4  internal 
waves  form  a  mode  3  internal  wave  with  the  frequency  (U3  and 
wave  number  ^3.  (b)  Equivalent  approach  as  in  Figure  4a  but 
applied  to  a  band  of  internal  wave  dispersion  relations  given  by  a 
range  of  geographical  latitudes  atxJ  vertical  stratification  profiles. 
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structure  at  higher  wave  numbers  differs  between  the  two  simulations  with  the  peak  shifting  from 
approximately  50  km  in  HYCOM12  to  35  km  in  HYCOM25. 

6.  Triad  Interaction  in  Wave  Number-Frequency  Space 

An  interpretation  for  the  structures  seen  in  E{K,  w)  (Figure  2)  and  UK,  to)  (Figure  3)  is  obtained  by  considering 
wave  triad  interactions  [McComas  and  Bretherton,  1 977],  where  we  assume  that  a  mode  1  internal  wave  with 
frequency  toi  and  wave  number  k,  interacts  with  a  mode  4  wave  (Figure  4a).  To  visualize  possibie  triad 
interactions,  the  mode  4  dispersion  relation  is  translated  in  K-co  space  by  as  denoted  by  the  green 
dashed  line.  Permitted  interactions  are  given  by  intersections  of  the  translated  dispersion  relation  with  the 
dispersion  relation  of  internal  waves.  As  an  example,  we  choose  the  intersection  with  the  mode  3 
dispersion  relation,  from  which  {k2,o)^  and  (k3/«3),  the  frequencies  and  wave  numbers  of  the  interacting 
and  resulting  internal  waves,  respectively,  are  then  specified.  By  definition,  the  resulting  wave  fulfills 
k3  =  k,  ±^2  and  <03  =  cu,  ±(i>2.  Note  that  Figure  4a  is  schematic  in  nature — for  simplicity,  we  have  not 
accounted  for  all  low-mode  interactions  nor  the  additional  complication  that  horizontal  wave  number  is  a 
vector  and  not  a  scalar  quantity. 

Interpretation  of  triad  interactions  in  K-w  space  can  be  extended  to  dispersion  relation  bands  instead  of 
lines  (Figure  4b).  The  bands  appear  because  of  the  range  of  f  and  N  in  the  considered  region  as 
described  in  the  previous  section.  The  intersections  consist  of  regions  which  appear  through  overlapping 
of  the  translated  dispersion  relation.  The  intersecting  regions  appear  as  the  patchy  structures  in  plots  of 
E{K,  co)  (Figure  2)  and  T(K,  ro)  (Figure  3).  Integration  over  these  fields  yields  a  wavey  rather  than  a 
monotonically  decreasing  E(co)  spectrum  (Figure  2)  especially  apparent  in  HYCOM12.  This  can  be 
interpreted  as  the  results  of  triad  interactions  of  a  limited  number  of  low-mode  internal  waves.  Since  our 
models  differ  in  horizontal  resolution,  patches  in  the  HYCOM12  and  HYCOM25  plots  differ  in  detail  but 
arise  from  the  same  physical  mechanism.  We  expect  that  the  strength  and  number  of  triad  interactions 
will  increase  with  finer  model  resolution. 

7.  Summary  and  Conclusion 

In  this  paper,  we  spectrally  analyzed  the  surface  horizontal  velocities  of  two  HYCOM  simulations  with  0.08° 
and  0.04°  horizontal  resolution,  respectively,  that  are  forced  simultaneously  by  atmospheric  fields  and 
astronomical  tides.  We  analyzed  a  two-dimensional  region  in  the  northeast  Pacific  north  of  the  Hawaiian 
Islands,  and  horizontal  wave  number-frequency  K-co  spectra  of  horizontal  kinetic  energy  and  nonlinear 
energy  transfers  are  computed  and  interpreted.  The  K-m  spectra  of  E  and  T show  clear  signatures  along  the 
dispersion  relations  of  the  low-vertical-mode  internal  waves.  The  internal  tides  are  visible  as  stripes  in  the 
K-co  spectra  and  peak  at  the  wavelengths  consistent  with  estimates  from  linear  internal  wave  theory. 
HYCOM25  has  larger  values  of  E  and  T  at  smaller  spatial  and  temporal  scales  than  does  HYCOM  12. 
Comparisons  of  model  frequency  spectra  with  seven  moored  current  meter  observations  suggest  that  the 
finer-resolution  HYCOM25  simulation  is  more  realistic  than  HYCOM12. 

We  conclude  from  this  study  that  calculating  the  spectral  nonlinear  kinetic  energy  transfer  in  K-co  space  is  a 
useful  technique  to  analyze  and  interpret  simulated  internal  wave  dynamics.  The  models  reproduce  low- 
vertical-mode  internal  wave  triad  interactions.  The  spectral  transfers  demonstrate  that  energy  is  extracted 
from  low-wave  number  near-inertial  and  supermode  3  tidal  rrtotions  in  order  to  feed  a  supermode  3  near- 
inertial  continuum  and  mode  2  2f  peak  in  the  internal  wave  spectrum.  The  patchy  structure  of  the  spectra 
can  be  interpreted  as  triad  interactions  of  a  limited  number  of  low-mode  internal  waves.  It  will  be  interesting 
to  follow  up  this  study  with  analyses  of  three-dimensional  model  fields  and  of  higher-resolution  simulations, 
which  are  likely  to  further  expand  the  realism  of  the  internal  wave  spectrum  in  global  ocean  models. 
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